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ABSTRACT 



Recent blazar detections by HESS, MAGIC, and VERITAS suggest that very- 
high-energy (VHE, E > 100 GeV) 7-rays may be produced in most, if not all, 
types of blazars, including those that possess intense circumnuclear radiation 
f-n* fields. In this paper, we investigate the interaction of nuclear VHE 7-rays with 

Q-r the circumnuclear radiation fields through 77 absorption and pair production, 

and the subsequent Compton-supported pair cascades. We have developed a 
Monte-Carlo code to follow the spatial development of the cascade in full 3- 
dimensional geometry, and calculate the radiative output due to the cascade 
as a function of viewing angle with respect to the primary VHE 7-ray beam 
(presumably the jet axis of the blazar). We show that even for relatively weak 
magnetic fields, the cascades can be efficiently isotropized, leading to substantial 



off-axis cascade emission peaking in the Fermi energy range at detectable levels 
for nearby radio galaxies. We demonstrate that this scenario can explain the 

q ■ Fermi flux and spectrum of the radio galaxy NGC 1275. 

O' 

Subject headings: galaxies: active -- gamma-rays: theory -- radiation mecha- 
nisms: non-thermal 



1. Introduction 

Blazars, a class of active galactic nuclei (AGNs) comprised of Flat-Spectrum Radio 
Quasars (FSRQs) and BL Lac objects, exhibit some of the most violent high-energy phenom- 
ena observed in AGNs to date. Their spectral energy distributions (SEDs) are characterized 
by non-thermal continuum spectra with a broad low-frequency component in the radio - 
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UV or X-ray frequency range and a high-frequency component from X-rays to 7-rays. They 
show rapid variability across the electromagnetic spectrum. In extreme cases, the very-high- 
energy (VHE ) 7-ray emission of blazars has been obser ved to vary on time scales of just a 
few minutes ( I Albert et al.l l2007al lAharonian et al.l 120071 ) . 



Leptonic and hadronic models (for a recent review see, e.g. lBottcherll2007l ) can generally 
successfully account for the overall SEDs observed from the known VHE 7-ray blazars, which 
are almost all BL Lac objects. BL Lac objects span a wide range of synchrotron peak 
frequencies, from IR to X-rays. According to the location of the synchrotron peak, they are 
classified as Low frequency peaked BL Lacs (LBLs: synchrotron peak in the IR), intermediate 
BL Lac objects (IBLs: synchrotron peak in the optical/UV) or High frequency peaked BL 
Lacs (HBLs: synchrotron peak in the X-rays). Until very recently, all known VHE 7-ray 
blazars were HBLs. However, the recent detections of the IB Ls W Comae, 3C 66A, and 
PKS 1424+240 by VERIT AS jAcciari et al.ll2008l bood. boioh. the LBLs BL Lacertae and 



S5 0716+714 by MAGIC (lAlbert et al. 



2007bl: lAnderhub et al-lboogh. an d even the Flat 



Spectrum Rad io Quasars (FSRQs) 3C 279 by MAGIC (lAlbert et al. 



20081 ) and PKS 1510- 



089 by HESS ( Wagner fc Beherall2010l ) suggest that most blazars are intrinsically emitters 
of VHE 7-rays. 

The fact that most LBLs and FSRQs have not been detected as VHE 7-ray sources 
might be primarily a result of the absorption of VHE 7-rays by lower-frequency (IR - optical 
- UV) radiation. As those objects tend to be located at greater cosmological distances 
(z > 0.2), 77 absorption on the Extragalactic Background Light ( EBL) becomes substantial 



suppressing any 



intrinsically produced > 100 GeV emission (e.g 



Dwek fc Krennrich 



Stecker fc Scully! l2008t iFranceschini et al.ll2008t iFinke fc Razzaqudl2009t iFinke et al. 



2005 



2010). 



However, possibly even more importantly, multi-GeV 7-rays produced in the high-radiation- 
density environment within the broad line region (BLR) and the dust torus of a quasar, is 



expected to be strongly at t enuated by 77 pair production 



Donea fc Protheroel 120031 : |Reimeij 12Q07J iLiu et al.l 12008 



e^ 



Protheroe fc Biermann 1997 



Sitarek fc Bednarek 



2008). The 



studies cited above show that for typical parameters expected in 3C279, photons above 
~ 100 GeV are expected to be strongly attenuated unless the 7-ray emission region is 
located beyond the broad line region. If VHE 7-rays are produced very close to the central 
engine, also 77 absorption by the direct accretion disk emission may become substantial 
rtSitarek fc Bednarek! Eoiol ). 



Given the now substantial number of non-HBL VHE 7-ray blazars, it seems plausible 
that most LBLs and FSRQs and their misaligned parent population, radio galaxies, are 
producing VHE 7-ray emission within their blazar zone. This view is also supported by 
the detection of VHE 7-ray emission from two non-blazar AGNs, namely the radio galaxies 
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M87 (jAharonian et al.ll2004j ) and Cen A (jAharonian et al.ll2009l ). The non-detection of VHE 
7-rays from most LBLs and FSRQs might then be due to the combined effects of local and 
intergalactic 77 absorption. As a consequence of 77 absorption in the local radiation field of 
the AGN, GeV - TeV electron-positron pairs are injected into the AGN environment, which, 
in the dense radiation field within the BLR and the dust torus of quasars, may initiate 
inverse-Compton supported pair cascades. 

The development of pair cascades induced by VHE 7-ray emission from blazars has so 
far concentrated on the development of Mpc-scale pair halos resul ting from the interactio n of 



(> 100 TeV) 7-rays with the Cosmic Micr owave backgro und (e.g. lAharonian et al.lll994[ ). or 
of VHE (> 100 GeV) 7-rays by the EBL (jVentersll2010l ). Due to the long Compton cooling 
timescale of pairs on the CMB or the EBL, those authors could reasonably consider the 
produced cascades isotropic, leading to an extended pair halo around the AGN. Small-angle 
deflection of seco ndaries in we ak in tergalactic magne tic fields has been included in those 
considerations by iPlagal (119951 ) and lElviy et al.l (120091). T h e dev elopm ent of pair cascades 
within the AGN has been discussed by iBednarek fc Kirk! (119951 ) and ISitarek fc Bednarek 

fcoioh . 



Depending on the magnetic field strength and orientation in the extended nuclear region, 
cascades developing within the high-radiation-energy-density environment within the BLR 
of a quasar may be efficiently isotropized in the immediate vicinity of the AGN. This will 
lead to distinct spectral features, which we will consider in this paper. In Section [2] we 
present some general considerations and analytic estimates of the expected results, including 
an estimate of the required magnetic field strengths within ~ a few pc from the central engine 
of a blazar (or radio galaxy), for which secondary electrons and positrons in pair cascades 
initiated by primary VHE 7-rays may be efficiently isotropized within the central region. 
The consequence would be quasi-isotropic radiation signatures from Compton-supported pair 
cascades, potentially observable in the Fermi energy range. In Section [3] we will describe a 
numerical code that treats the full three-dimensional development of these cascades, together 
with the general model setup and simplifying assumptions. Numerical results for generic 
parameters will be presented in Section HJ In Section [51 we will d emonstrate that the recent 
Fermi detection of the radio galaxy NGC 1275 (jAbdo et al.ll2009a ) can be plausibly explained 
by the cascade emission from a misaligned VHE 7-ray emitting blazar. We summarize and 
present an outlook towards future work in Section |6j 
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2. General Considerations and Estimates 

In the limit of a VHE 7-ray photon with energy E 1 interacting with an IR/optical/UV 
photon from the dust torus or the BLR with energy E s <C E^, an electron-positron pair 
with particle energy E e = 7 m e c 2 ~ E 1 /2 = 1 E Te y TeV will be produced, moving along the 
direction of the primary 7-ray photon to within an accuracy of A9 ~ I/7. Depending on the 
strength of the magnetic field near the point of pair production, B = 1 B-q fiG, they will be 
deflected on a length scale of the order of the Larmor radius, 

r g ~lQ- 3 ET*vB^pc. (1) 

To investigate whether particles are efficiently isotropized before producing secondary syn- 
chrotron and/or inverse-Compton (IC) emission, the isotropization length has to be com- 
pared with the radiative cooling length, A sy and Aic, respectively The synchrotron cooling 
length can be estimated to 

A sy ~ 3.8 x 10 6 E^ v BZI pc. (2) 

To estimate the IC cooling length on an external radiation field from the BLR we calculate its 
energy density as wblr ~ -^D r BLR/(4vr_R| LR c), which we parameterize through the accretion 
disk luminosity L D = 10 46 L 46 erg s _1 and a BLR with a reprocessing optical depth of 
tblr = 0.1 r_i at an average distance of -Rblr = 0.1 R-i pc from the central engine. Hence, 
the total luminosity of the BLR will be Lblr = TblrAd- If Compton scattering occurs in 
the Thomson regime, we find 

Aic ~ 5 x 1(T 6 E£ v L^tzIR\ pc. (3) 

This illustrates that one may expect particles at energies substantially below 1 TeV to be fully 
isotropized, while higher-energy particles might lose a substantial fraction of their energy 
while still traveling along the primary VHE 7-ray beam. Comparison of the synchrotron and 
IC cooling lengths suggests that the radiative output from the se condaries will be strongly 



domi nated by IC emission, initiating an IC-supported cascade (e.g. |Protherodll986l ; IZdziarski 



19881 ). At low frequencies far below the pair production threshold, produced by secondaries 
which are fully isotropized (i.e., Aic ^ r g)-, the cascade spectrum will obtain a vF v oc v 1 ! 2 
shape. This is the consequence of the secondaries being injected at high energies and then 
being subject to Compton cooling in the Thomson regime, resulting in a -AT (7) oc 7~ 2 pair 
spectrum. This low-energy spectral shape (though not its total flux) will be independent 
of the primary gamma-ray spectrum. Two effects will produce a turnover towards higher 
frequencies. First, for any given viewing angle 9 with respect to the direction of propagation 
of the primary 7-ray (being absorbed in the 77 pair production process), we can find a critical 
electron energy for which the deflection angle over a Compton length equals the observing 
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angle, i.e., 9 ~ \\c/ r g- Higher-energy particles will radiate preferentially at smaller viewing 
angles, while lower-energy particles can efficiently contribute to the emission at the given 
angle. This yields the characteristic electron energy E e ^ T corresponding to a given observing 
angle 6: 



E eM = m e c 2 J - 3eB - ~ 70 B% 2 A.! L^ 2 rZl l% IT 1 ' 2 GeV. (4) 

V 4(7 T M B LRP 

If these electrons can scatter the soft photon field in the Thomson regime, a turnover 
should occur at photon energies of 

E lCM = . 3eB , E s ~ 18 5_ 6 R\ Lll tZ\ 9- 1 (%) GeV. (5) 

We point out that, in this estimate, we assumed that the Compton cooling length is a 
realistic measure of the distance travelled by the electron/positron since its production from 
annihilation of a primary 7-ray. This approximation will break down if the particles require 
multiple Compton scatterings to reach the break energy corresponding to Eq. [5j This is 
due to both an increasing Compton cooling length and a decreasing Larmor radius as the 
electron energy is reduced by previous scatterings. 

Second, if Compton scattering to energies given by Eq. [5] occurs in the Klein-Nishina 
regime (which is the case if £ic,br > ^e,br) a turnover is expected at the transition from 
Thomson to Klein-Nishina scattering at 

£ic,kn = 260 (^\ GeV. (6) 



3. Model Setup and Code Description 

The geometrical setup of our model system is illustrated in Figure HJ The primary 
VHE 7-ray emission from the blazar zone is represented as a mono-directional beam of 7- 
rays propagating along the X axis. The incident VHE 7-ray spectrum is represented by a 
straight power-law with photon index a. Those 7-rays may interact via 77 absorption and 
pair production with a radiation field. For a first general investigation and feasibility study 
presented here, we approximate the radiation field as monoenergetic and isotropic within a 
fixed boundary, given by a radius -R ex t, i-e., 
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Observer 




Fig. 1. — Geometry of the model setup. 



w ext (e, r, Q) = u 8(e - e cxt ) H(R ext - r) (7) 

where H is the Heaviside function, H (x) = 1 if x > and if (x) = otherwise. A magnetic 
field of order ~ /zG is present. Without loss of generality, we choose the y and z axes of our 
coordinate system such that the magnetic field lies in the (x,y) plane. 

We have developed a Monte-Carlo code which treats the processes of 77 absorption and 
pair production, 7-ray and electron/positron propagation, and Compton scattering. The 
code generates a single 7-ray photon at a time, at the origin of our coordinate system, 
propagating in the x direction. In order to improve the statistics of the otherwise very few 
highest-energy photons, we introduce a statistical weight inversely proportional to the pho- 
ton energy. The code calculates the absorption opacity k 77 , using the full pair production 
cross section. Based on the corresponding absorption length, a location for the occurrence 
of the next 77 absorption / pair production process is drawn. If the pair production site 
is outside the radius -R ex t, the photon escapes; otherwise, the photon is absorbed, and an 
electron-positron pair is created. The code uses the analytic result for the 77 pair produc- 



tion spectrum of iBottcher fc Schlickeiserl (119971 ) to draw the energies of the electron and 
the positron. As motivated above, we assume that both particles initially propagate in the 
direction of the absorbed 7-ray For both particles, the Compton scattering length Aic is 
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calculated using the full Compton cross section. As we expect the magnetic-field energy 
density to be much smaller than the radiation energy density, we neglect synchrotron losses 
to the electrons/positrons. Based on the value of Aic, the code draws a length that the 
electron/positron travels before the next Compton scattering event occurs. It then propa- 
gates the electron/positron, using the full 3-D geometry, through its gyrational motion in 
the magnetic field, to calculate the (x,y,z) coordinates and direction of motion of the elec- 
tron/positron at the point of scattering. If this point of scattering is outside the radius -R ex t, 
the electron/positron escapes; otherwise, Compton scattering occurs. The energy of the pro- 
duced photon is drawn us ing a 5(Q SC — Q e ) approximation for the Compton cross section (e.g., 



Dermer fc Bottcherll2006l ). i.e., the scattered photon is traveling in the same direction as the 
electron/positron before scattering. The produced (7-ray) photon is then tracked through 
the same photon tracking routine as the primary VHE 7-ray photons (properly accounting 
for the correct location and direction of propagation). The energy of the electron/positron 
is reduced by the energy of the scattered photon, and the particle is returned to the pair 
tracking routine. If the energy of the electron/positron is below a set threshold (determined 
by the condition that they will no longer produce Compton emission in the energy range of 
interest), the code will move on to the next particle. 

The energies, statistical weights, and directions of propagation of photons escaping from 
the region of high external radiation field (bounded by R ex t), are written into a photon event 
file. In a post-processing routine, this event file will be read to produce photon spectra with 
arbitrary energy and angular binning. 



4. Numerical Results 

We have used the cascade Monte-Carlo code described in the previous section to evaluate 
the angle-dependent cascade spectra for a variety of generic parameter choices. Figure [2] 
illustrates the viewing angle dependence of the cascade emission. For this simulation, we 
assumed a magnetic field of B = 1 /iG, oriented at an angle 9b = 5.7° with respect to the X 
axis (B x = 1 /zG, By = 0.1 /iG). The external radiation energy density is uq = 10 -2 erg cm -3 , 
extended over a region of radius i? ex t = 10 19 cm with photon energy E s = E^ ya . The 
incident 7-ray spectrum has a photon index of a = 2.5 and extends out to -E^max = 5 TeV. 
The results have been normalized to a flux level in the forward direction corresponding to 
a 7-ray bright blazar. The spectra for all other directions have been normalized with the 
same normalization factor. The curves are labeled by the cosine of the observing angle, 
/1 = cos^obs- In the forward (blazar) direction, one clearly sees the 77 absorption cut-off at 
an energy E c = (m e c 2 ) 2 /E s ~ 25 GeV. The cutoff is very sharp in this simulation because of 
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Fig. 2. — Cascade emission at different viewing angles (/i = cos# bs)- Parameters: B = 1 /zG, 
9 B = 5.7°; u = lO" 2 erg cm" 3 , R cxt = 10 19 cm, E s = E Lya , a = 2.5, £ 7 , max = 5 TeV. The 
solid (red) straight line indicates the incident primary 7-ray spectrum. 



our S approximation of the external radiation field, combined with a high 77 absorption depth 
near threshold for the parameters chosen here. Below this cutoff, the forward-component of 
the cascade emission leads to a slight bump beyond the primary 7-ray power-law spectrum. 
The cascade emission at larger viewing angles \i < 0.85 has a low-frequency shape close to 
the expected vF v oc z/ 1//2 behaviour, and exhibits the progressive suppression of the cascade 
emission at high energies with increasing viewing angle due to incomplete isotropization of 
the secondary particles at high energies, as expected from Eq. [51 

Figure |3] illustrates the effect of a varying external radiation field energy density uq. For 
energy densities u > 1CT 3 erg cm -3 , 77 absorption is essentially saturated, i.e., all VHE 
photons above the pair production threshold will be absorbed. Hence, the magnitude of 
the cascade becomes almost independent of uq. For smaller w , the decreasing flux in the 
cascade emission reflects the decreasing fraction of VHE 7-ray photons absorbed. This latter 
is the regime in which our estimate of the turnover frequency (Eq. |5]) is applicable. In the 
saturated regime, the electrons/positrons have to undergo many scatterings before reaching 
the isotropization energy so that the Compton scattering length is no longer an appropriate 
measure of the distance traveled, as assumed in the derivation of Eq. |5j As expected, the 
low mq case results in a larger turnover energy than the high uq cases. 
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Fig. 3. — The effect of a varying external radiation energy density Parameters are the same 
as for Figure |2] in the angular bin 0.38 < \i < 0.54. 



The effect of a varying magnetic field strengh for fixed magnetic field orientation (9b = 
45°) is illustrated in Figure SJ We see that the cascade development is extremely sensitive 
to the transverse magnetic field B y for weak magnetic fields. The cascades exhibit a rapid 
transition to the limit in which even the highest-energy secondary particles are effectively 
isotropized before undergoing the first Compton scattering interaction. Hence, for magnetic 
field values expected on sub-pc or pc scales around an AGN (B ~^> 1 nG) and large angles 
9b, there is no pronounced break in the cascade spectrum out to large energies near the 77 
absorption trough at E c . 

Figure |5] shows the effect of a varying magnetic-field orientation with respect to the jet 
axis, for fixed magnetic- field strength B = 1 /xG. The figure illustrates that it is primarily the 
perpendicular (B y ) component of the magnetic field which is responsible for the isotropization 
of secondaries in the cascade. Obviously, for a perfectly aligned magnetic field (9b = 0°), 
our code does not predict any cascade emission in any off-axis direction, since we neglect 
the spreading of the cascade due to the kinematics of the pair-production process or recoil 
from Compton scattering. For a small inclination angle, the small perpendicular magnetic- 
field component leads to inefficient isotropization and, hence, a break in the off-axis cascade 
spectrum at low energies. For large inclination angles (B y > B x ), isotropization becomes 
very efficient out to energies close to the 77 absorption trough. 
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Fig. 4. — The effect of a varying magnetic field strength, for a fixed angle of 9b = 45° 
between jet axis and magnetic field. All other parameters are the same as for Figure |2] in 
the angular bin 0.38 < [i < 0.54. 



5. Application to Radio galaxies: The case of NGC 1275 



From the numerical results shown above it has become obvious that VHE 7-ray induced 
cascades even in the presence of rather weak (~ /iG) magnetic fields can be efficiently 
isotropized and produce a substantial MeV - GeV 7-ray flux in directions misaligned with 
respect to the jet axis (the forw ard or "blazar" direction). The standard AGN unification 
scheme (jUrry fc Padovanil Il995l ) proposes that blazars and radio galaxies are intrinsically 
identical objects viewed at different angles with respect to the jet axis. According to this 
scheme, FR I and FR II radio galaxies are believed to be the parent population of BL Lac 
objects and FSRQs, respectively. Hence, if most blazars, including LBLs and FSRQs, are 
intrinsically VHE 7-ray emitters potentially producing pair cascades in their immediate 
environments, the radiative signatures of these cascades might be observab l e in ma ny radio 
galaxies. I n fact, already three radio galax i es (NG C 1275: lAbdo et al.l (12009al ). M 87: 



Abdo et al.l (j2009bl ). and Cen A: lAbdo et al.l ( ]2009d )) have been detected by Fermi, while 
EGRE T provided evidence for > 100 MeV 7-ra y emission from two m ore radio galaxies, 
3C 111 ^Nandikotkur et al-lbooA and NGC 6251 JMukheriee et al.l l20o3). In this paper, we 
focus on the radio galaxy NGC 1275 ( lAbdo et al.ll2009al ) and investigate whether the Fermi 
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Fig. 5. — The effect of a varying magnetic field orientation, for a fixed magnetic field strength 
of B = 1/xG. All other parameters are the same as for Figure [2] in the angular bin 0.38 < 
fj, < 0.54. 



spectrum and flux level are consistent with the interpretation of cascade emission from a 
misaligned, 77 absorbed blazar. 

NGC 1275 is located at a distance of d = 74 Mpc. It hosts the powerful FR I radio galaxy 
Perseus A (= 3C 8 4). The radio jet is directed at an angle of 9 « 30° - 55° with respect 
to the line of sight ([Walker et al.lll994f ). NGC 127 5 also hosts a Seyfert -like nucleus with a 
total line luminosity of Lblr = 1-6 x 10 42 erg s _1 (jZirbel &: Baumlll995l ). Fermi observed a 
7-ray flux of F(> 100 MeV) = (2.10 ±0.23) x 10~ 8 ph cm" 2 s" 1 from NGC 1275, w hich could 
be fit by a power-law with a photon index of T — 2.17 ± 0.05 (jAbdo et al.ll2009al ). However, 
the Fermi spectrum does seem to exhibit a marked steepening towards higher energies. A 
comparison with the EGRET upper limit indicates substantial flux variability on time scales 
of years to decades, indicating that the 7-ray emission must be produced in a region of R < a 
few pc. 

Our simulation results can be properly normalized to the expected flux for a radio galaxy 
with the following considerations: The forward spectrum shown in Figure [2] corresponds 
to the escaping radiation in a solid angle interval AQ = In A/x, where A/j = 0.15 in the 
examples shown here. However, the cone of blazar emission might be smaller, namely Afi& w 
7T# 2 , where 9b is the opening angle of the blazar emission cone. Hence, a realistic flux 
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Fig. 6. — Fit to the Fermi spectrum of NGC 1275 with a simulated cascade spectrum from 
a mis-aligned blazar, along with the cascade spectra at larger viewing angles. 



normalization for any off- axis direction requires a correction factor of / ang = 0%/(2 A/i). The 
normalization will also depend on the difference in distances, /dist = (d h \/d RG ) 2 , where du 
and g?rg are the luminosity distances of the blazar and the radio galaxy, respectively. Hence, 
we apply a total normalization factor / = / ang /dist- 

Figure |6] illustrates that the Fermi spectrum of NGC 1275 can be well matched with a 
VHE 7-ray induced cascade from a misaligned blazar. For the normalization of the incident 
7-ray beam, we used the Fermi spectrum of the blazar 3C279 at dj, = 2.2 Gpc. The external 
radiation field is parameterized through w ext = 5 x 10~ 2 erg cm -3 and R r ^ t = 10 16 cm. This 



size s cale is appropriate for low-luminosity AGN as observed in NGC 1275 (e.g. Kaspi et al. 
20071 ). and the parameters combine to a BLR luminosity of Lblr = 47ri?L+ c ?w = 1.9 x 



''ext 



10 42 erg s _1 , in agreement with the observed value for NGC 1275. The magnetic field is 
B = 1 mG, oriented at an angle of 9b = 8°. The cascade spectrum shown in Figure [6] 
pertains to the angular bin 0.73 < \i < 0.87 (corresponding to 30° < 6 < 43°), appropriate 
for the known orientation of NGC 1275. The factor required to re-normalize between the 
simulated spectrum in the forward direction (at a flux level observed for 3C 279) and the 
sideways spectrum, fit to NGC 1275, was / = / ang /dist = 146. This corresponds to a blazar 
cone opening angle of 9\, = 12.7° for 3C279. It should be pointed out that our normalization 
to flux levels of 3C279 are only meant to demonstrate that our parameter values correspond 
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to conditions reasonably expected in misaligned blazars in general. The numbers quoted 
above pertaining to 3C 279 should not be taken at face value as a diagnostic for 3C279 itself. 



6. Summary and Outlook 

We have developed a Monte-Carlo code to simulate the full 3-D development of VHE 
7-ray induced pair cascades in the weakly magnetized environments of blazars and radio 
galaxies. These cascades develop as Compton-supported pair cascades if there is a substantial 
radiation energy density, e.g., from the broad-line region or the dust torus around the central 
engine. We have demonstrated the angle-dependence of these cascades as a function of the 
magnetic field and the radiation energy density in the AGN environment. In the off- axis 
directions, these cascades exhibit a vF v oc v x l 2 spectrum at low energies and a turnover 
towards higher energies due to incomplete isotropization of high-energy secondares in the 
cascade, depending on the strength and orientation of the magnetic field and the radiation 
energy density. 

We have demonstrated that the off-axis cascade emission may well be detectable from 
nearby radio galaxies, which, according to the blazar unification scheme, are the mis-aligned 
parent population of blazars. We have presented a fit to the Fermi spectrum of the radio 
galaxy NGC 1275 with parameters expected if this radio galaxy is the mis- aligned counterpart 
of a typical Fermi-detected blazar. 

The present work is to be understood as a general proof of the potential importance of 
VHE 7-ray induced pair cascades in blazar environments. For this purpose, we have used a 
simplified representation of the radiation field as monoenergetic, homogeneous and isotropic 
out to a limiting radius. In future work, we will relax this assumption to allow for an 
arbitrary radiation spectrum and more realistic spatial and angular dependence. While we 
expect slight quantitative differences in the spectral shape of the emanating cascade emission 
due to these modifications, the general viability of the process as well as its dependence on 
magnetic field properties and the external radiation field are expected to remain robust 
predictions of our work. 

We thank Jun Kataoka for sending us the Fermi data points for NGC 1275. This work 
was supported by NASA through Fermi Guest Investigator Grant NNX09AT81G. 
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